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Abstract 

We study the Thirring model in three spacetime dimensions, by means of Monte 
^ ' Carlo simulation on lattice sizes 8"^ and 12^, for numbers of fermion flavors Nf — 

2, 4, 6. For sufficiently strong interaction strength, we find that spontaneous chiral 



O 

^ symmetry breaking occurs for Nf = 2,4, in accordance with the predictions of 

^ \ the Schwinger-Dyson approach. The phase transitions which occur are continuous 

. and with critical scaling behaviour depending on Nf. For Nf = 6 our results are 

p preliminary, and no firm conclusions about the existence or otherwise of chiral 
CD ! 

r-| , symmetry breaking are possible. 
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The three dimensional Thirring model is a field theory of relativistic fermions interact- 
ing via a contact term between conserved vector currents. Its Lagrangian is written 

C = + 171)11;, + ^(V'^7M^^)^ (1) 

where ^ire four-component spinors, m is a bare, parity-conserving mass, and the 

index i runs over Nf distinct fermion species. Since the coupling has mass dimension 
— 1, naive power-counting suggests that the model (1) is non-renormalisable. However, as 
has been suspected for many years [1,2], an expansion in powers of 1/A^/, rather than g'^, 
is exactly renormalisable. At leading order in l/Nf, interaction between vector currents 
is dominated by exchange of a fermion - anti-fermion bound state described in terms of 
a chain of vacuum-polarisation "bubbles" . In the ultra-violet limit the interaction is thus 
transformed from a monmentum-independent contact term to a softer A/ {Nfk) behaviour, 
where ^ is a numerical constant independent of g: this asymptotic behaviour can be used 
to evaluate divergent graphs at higher orders in l/Nf, eg. in [3], where renormalisability 
of the massless model is explicitly demonstrated to 0{1/Nf). 

The property of renormalisability signals that the model's expansion exhibits a 

UV-stable fixed point of the renormalisation group, the continuum limit being taken in the 
limit g^A oo, where A is a UV cutoff. RG fixed points have also been observed in other 
three-dimensional four-fermi models [4,5] . The distinctive feature of the Thirring model is 
that for d < 4 the vacuum polarisation is UV-finite so long as the regularisation respects 
current conservation (eg. Pauli- Villars) . This means that there is no need to fine-tune 
g'^ to a critical value: a continuum limit may be taken for any value of the dimensionless 
parameter mg^ (at least to leading order in [3]), the theory thus obtained having a 

variable ratio of, say, physical fermion mass to vector bound state mass. In the RG sense 
the interaction {'i/jji^if;)'^ is a marginal operator for 2 < d < 4, whereas, say, the interaction 
{iptpY in the Gross- Neveu model is relevant [6]. 

Another possibility raised by the expansion is the equivalence of the Thirring 
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model in the strong-coupling limit mg^ — > oo, in which the vector particle becomes mass- 
less, with the infra-red limiting behaviour of QED in three spacetime dimensions. In 
massless large- A^/ QED3 vacuum polarisation screens one-photon exchange to the extent 
that the interaction is again transformed to 1/k [7]. The 0(1 /Nf) corrections to 

the models evaluated in respectively UV (Thirring) or IR (QED) limits appear to coincide 
[3,8]. 

The 1/Nf expansion may not, however, describe the true behaviour of the model, 
particularly for small Nf. For instance, spontaneous chiral symmetry breaking, signalled 
by a vacuum condensate (V'V')) is forbidden to all orders in 1/Nf, and yet may be predicted 
by a self-consistent approach such as solution of the Schwinger-Dyson equations [2,9,10]. In 
this approach a non-trivial solution for the dressed fermion propagator S{p) = {A{p)i;^ + 
E(p))~-^ is sought, ie. one in which the self-energy E(p) and hence {'ipip) are non-vanishing 
in the chiral limit m ^ 0. Unfortunately, the SD equations can only be solved by truncating 
them in a somewhat arbitrary fashion. The usual approximation [2,10] is to assume that 
the vector propagator is given by its leading-order form for m = in the 1/A^/ expansion, 
viz. 

and that the fermion- vector vertex function is well-approximated by the bare vertex (the 
so-called "planar" or "ladder" approximation): 

r,{p,q) = -^i,. (3) 

The longitudinal part of D/^l, raises potential ambiguities: the most systematic treat- 
ment has been given by Itoh et al [10], who note the equivalence of the Thirring model 
with a gauge-fixed form of a fermion-scalar model possessing a local gauge symmetry and 
then use a non-local gauge-fixing condition to find a gauge in which the "wavefunction 
renormalisation" A{p) = 1. This simplification enables the SD equations to be exactly 
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solved in the limit — > oo, with the result that a non-trivial solution for exists for 

Nf<Nfc=^c^ 4.32. (4) 

Moreover since the integral equations require the introduction of a UV cutoff A, a feature 
of this solution is that the induced physical mass scale depends on Nf in an essentially 
singular way: 

oc exp I , ^ I : ocA^a^ oc exp ( , \ . (5) 

This implies that a continuum limit only exists as Nf Nfc, the scenario being very 
similar to that proposed by Miranskii and collaborators for strongly-coupled QED4 [11]. 
Unfortunately no analytic solution exists for (7^ < 00; however using different techniques 
Kondo [12] has argued that a critical line Nfc{g'^) exists in the {g'^,Nf) plane, which is a 
smooth invertible function. Therefore for integer Nf < Nfc one might expect a critical 
scaling behaviour 



/ a 



,2 



\ vfe-1 



(6) 



corresponding to a symmetry restoring transition at some critical point g = g^. Pre- 
sumably in this scenario the Thirring interaction has become relevant: there may exist 
a novel strongly-coupled continuum limit at the critical point not described by the 1/A^/ 
expansion. 

There are good reasons to be cautious of this picture, however. Using a different 
sequence of truncations Hong and Park have found chiral symmetry breaking for all Nf 
[9], with 

1 / NfTT 



-^ocexp(^-^j, (7) 

a result which is non-analytic in 1/Nf. Moreover in the limit g^ ^ 00 the system of 
SD equations obtained are very similar to those of large- A^/ QED3, in which case studies 
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beyond the planar approximation, using improved ansatze for the vertex F^, suggest that 
the condition A{p) = 1 is unphysical, and that chiral symmetry is spontaneously broken 
for all Nf [13]. In the current context this would imply < oo for all Nf. 

For these reasons we consider a numerical study of the lattice-regularised model to 
be timely. If, as suggested above, the Thirring model lies in the same universality class 
as QED3, then a numerical study may shed light on the value of Nfc for this model 
[7]; previous lattice studies [14] have been plagued by large finite volume effects due to 
the slow fall-off of the photon propagator cx: 1/x. The corresponding propagator in the 
Thirring model falls as 1/x^, so sensible results may emerge on smaller systems. A second 
motivation is the possible existence of a novel continuum limit: since vacuum polarisation 
corrections to the vector propagator are finite to all orders in 1/A^/ there should be no 
competing effects of charge screening, which obscures the issue in QED4. 

The lattice action we have used is as follows: 

x^ii xi 

+ -^Yxi{x)x^{x + (i)Xj{x + lj)Xj{x) 
= \ Yxiix)v^^ix) {1 + iA^{x)) Xi{x + fi) + h.c. m ^ Xi(^)Xi(^) + ^^^^U^)^ 

x^i xi xfj, 

= ^Y1 Xi{x)M^A,m){x,y)xi{y) + 

xyf/,i xn 

(8) 

where x, X are staggered fermion fields, r]^ the Kawamoto-Smit phases, m is the bare mass, 
the flavor index i runs from 1 to N , and we have introduced for the fermionic 

bilinear, which depends on both the auxiliary field and the mass. The second form of 
the action is the one actually simulated: the equivalence of the two forms follows from 
Gaussian integration over the real- valued auxiliary field defined on the lattice links (for 
N = 1 there is an alternative formulation in terms of a compact complex-valued auxiliary 
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[15,16]). The vector-like interaction of the action allows the introduction of a checker- 
board, which in turn enables simulation of the system for any N. In three Euclidean 
dimensions staggered fermions describe two continuum species of four-component fermions, 
with a parity-conserving mass term [17]. Hence the number of physical flavors Nf = 2N. 
An interesting feature of the lattice formulation (8) is that the interaction current is not 
exactly conserved. The conserved current in lattice gauge theory incorporates the gauge 
connection exp(z^^). This means that at leading order in 1/N the vector propagator 
receives an extra contribution from vacuum polarisation, essentially due to the absence of 
the diagram of Fig. 1. The effect can be absorbed into a redefinition of the coupling: 

^^= i-,'jM ' 

where J(m) is the value of the integral depicted in Fig. 1. The physics described by 
continuum 1/Nf perturbation theory occurs for the range of couplings gj^ G [0, oo), ie. for 
e [Q,giim)'i leading order in 1/N 

= Jim); with J(0) = ^. (10) 

9lim ^ 

We therefore expect to see some kind of discontinuous behaviour in our simulations for 
small values of l/g^. 

In this letter, we aim to clarify the chiral symmetry breaking pattern by studying the 
chiral condensate: 

which, in the limit m ^ 0, is an order parameter for the spontaneous symmetry breaking. 
We performed simulations on 8^ and 12^ lattices for Nf = 2, 4, 6, using the hybrid Monte 
Carlo algorithm. Bare mass values m ranged from 0.4 down to 0.02, with most attention 
paid to the range 0.05 - 0.02. For each mass and coupling we performed roughly 500 HMC 
trajectories, the trajectory length being drawn from a Poisson distribution with mean 0.9. 
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The condensate {'ipip) was measured with a stochastic estimator every few trajectories. 
To maintain reasonable acceptance rates the timestep varied from 0.15 on 8^ at m = 0.4 
down to 0.022 on 12^ at m = 0.02. We found that considerably more work was needed 
to perform matrix inversion in this model than for the Gross-Neveu simulations described 
in [5]. Another difference is that in this case since the critical region of interest occurs 
at successively stronger couplings as Nf is raised, the CPU required also grows with Nf, 
despite the 1 /Nf suppression of quantum corrections. 

In Fig. 2, we plot a = {ipip) vs. l/g"^ for the three values of Nf studied, for m = 0.10 
on a 8^ lattice. The models with different Nf have apparently coincident condensates in 
the strong-coupling region l/g"^ < 0.3, but thereafter the (i^i^) signals peak to maxima at 
distinct values oil/ before falling away. It is tempting to associate the strong-coupling 
region with < from the discussion following Eq. (9), although the correspondence 
with the value of gf^^ predicted in Eq. (10) is not good. It may well be that the value of 
the diagram of Fig. 1 is considerably altered in a chirally broken vacuum. 

As stated above, in order to study spontaneous chiral symmetry breaking, one has to 
monitor the value of as m ^ 0. Our results for the chiral condensate for Nf = 2 

are reported in Fig. 3 for different values of the bare mass. A naive extrapolation to the 
chiral limit from the lattice data at fixed l/g"^ is probably unreliable in the range of bare 
masses we have explored. In order to determine the critical point, we need to perform a 
global fit of our data incorporating many values of m and l/g^. Therefore, we have to use 
an equation of state (EOS) relating the external symmetry breaking parameter m to the 
response of the system (i/^'iIj) ^tnd the coupling 1/g^ [16,18]. 

A generic EOS, inspired by the critical behaviour of spin systems, can be written in 
terms of the scaled variables: 

m(V'V)-' = ^(A(l/5')(V^)^/^), (11) 
where A(l/^^) = 1/gl — 1/g'^ is the reduced coupling. At g = gc, Eq. (11) is the usual 
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scaling relation: 

while a Taylor expansion for small A(l/(7^) yields: 

m = B{ijijf + AA(1^2)(^^)*-i//3 + . . . (12) 

where A, B = J-''{0),J-'{0) respectively. At this stage, one sees that, for vanishing m, 
Eq. (12) is simply the definition of the critical exponent /?: 

and that there are no logarithmic corrections, since these should only appear in 4-d [18]. 
If the critical behaviour is described by mean- field theory, we expect S — 3 and (3 — 1/2, 
yielding: 

(«' = «i7^ + «2A(l/^^), (13) 

which shows that {'ipip)'^ is a linear function of the ratio m/{'ip'ip). Such a plot is known 
as Fisher plot. From Eq. (13) we see that a positive value of the intercept corresponds 
to a non- vanishing value of the chiral condensate for m = 0, while the intercept will be 
exactly zero at the critical coupling. In Fig. 4 and 5 we show the Fisher plots for A^/ = 2, 4 
respectively, where we can see at first glance an indication of chiral symmetry breaking, 
according to the criterion stated above. In order to get a more quantitative evidence, we 
have fitted our data using Eq. (12) and a simpler version of it based on the hypothesis that 
5 — 1/P = 1 [19], which we will call respectively fit I and II in what follows. We should 
stress that in the absence of a systematic critical theory the forms I and II are used simply 
as effective descriptions of the data. The values of the fit parameters, their errors and the 
are listed in Table 1. The number of values of the chiral condensate included in the fit 
is chosen in order to minimize the value of the reduced x^. The results of fit II are shown 
in Figs. 4 and 5 and seem to describe the data quite well. The dashed line in Fig. 3 is the 
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curve one obtains using the results of fit II with m = 0. It shows clearly that, within the 
range of m we have explored, the chiral condensate is still far from its chiral limit value, 
thus providing a justification a posteriori for the impossibility of extrapolating the data 
naively. 

There are a few conclusions one can draw from the numerical analysis that we would 
like to emphasize. First, for both values of Nf, we find clear evidence of chiral symmetry 
breaking at finite values of the coupling, as predicted by the Schwinger-Dyson approach. 
From Eq. (7) we get: 

|^^exp,-.VS)=0...1 

which is not too far from the fitted value 0.342 ±0.015 (using the data from fit II). 

Secondly, although any claims to understand the details of the critical scaling must be 
premature, the fits strongly suggest that the models with Nf = 2 and 4 are described by 
distinct critical theories, in a sense on "either side" of the mean field theory. If we relax 
the requirement 5 — 1/ (3 = 1 (which means using fit I instead of fit II) then the difference 
in the fitted values of 5 becomes even more apparent. This is significant because similar 
EOS fits in QED4 reveal no such differences between Nf = 2 and Nf = A [20] . 

Finally, we report some preliminary results for Nf = 6. Figure 6 shows ('0V') vs. m 
for two values of g. The = 0.5 data suggest a linear extrapolation to the chiral limit, 
yielding a small condensate equal to 0.013(4). For l/g^ = 0.4 it is less clear how to make 
the extrapolation. Clearly in either case reliable data at much smaller mass values would 
be needed for confirmation or otherwise of chiral symmetry breaking: comparing data from 
different lattice sizes, we have found that finite size effects become more important as we 
go to larger Nf, which means that larger lattice sizes will probably be needed before we 
can proceed to a more quantitative study. 

We conclude by briefiy summarizing our results. We have shown by numerical simu- 
lations that spontaneous chiral symmetry breaking does occur in the Thirring model for 
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finite Nf, in contradiction with the perturbative expansion, but in agreement with 

the Schwinger-Dyson approach at least for Nf ~ 2,4. In these two cases we were able to 
determine the critical coupling and the critical exponent S by fitting to a plausible EOS, 
and the fits suggest that the two models have different critical behaviour. For Nf = 6, we 
were unable to find clear evidence in favour of symmetry breaking, but cannot yet exclude 
a non- vanishing condensate in the chiral limit. In the future we plan to investigate in more 
detail the theory for Nf = 2, 4, 6 at the critical point, focussing on critical exponents, the 
renormalized charge, and spectroscopy. 
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par am. 


fit I 


fit II 


Nf = 2 


Pc 


2.03(9) 


1.94(4) 


S 


2.32(23) 


2.68(16) 


P 


0.71(9) 




A 


0.32(5) 


0.37(1) 


B 


1.91(43) 


2.86(53) 


xVd.o.f 


2.4 


2.1 


Nf = A 


/3c 


0.63(1) 


0.66(1) 


5 


3.67(28) 


3.43(19) 


P 


0.38(4) 




A 


0.78(5) 


0.73(2) 


B 


7.9(2.8) 


6.4(1.5) 


xVd.o.f 


3.1 


2.0 



Table 1 

Results from the fits 
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Figure Captions 

Figure 1: Diagram contributing to coupling constant renormalisation at leading order in 

Figure 2: chiral condensate u vs. 1/^^ for N = 1,2,3, corresponding respectively to 
Nf = 2,4,6. 

Figure 3: chiral condensate a vs. for Nf = 2 and different values of the bare mass m. 

Figure 4: Fisher plot for Nf = 2, from data at /? = 1.6(A), 1.8(<), 2.0(V), 2.2(>), 2.4(+), 
2.6(x). 

Figure 5: Fisher plot for Nf = 4, from data at /? = 0.5(o), 0.6(.), 0.7{o), 0.8(A), 0.9(<), 
l.O(V), 1.1(>), 1.2(+), 1.3(x), 1.4(*). 

Figure 6: chiral condensate a vs. m for different values of the coupling (here P = 
and is not related to the critical exponent). 
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Figure 1 
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Fig. 2 
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